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AB8IBACT 


Guided  eleetroaagDetie  nodes  elogog  a  homogeneous  circular  columa  of 
plasna  lirit>edded  in  free  space  are  examined  for  the  two  extreme  eases  of 
sero  and  large  axial  nagnetie  field.  The  plasma  is  assuawd  to  he  cold 
and  coUlsioinless  and  loo  notion  is  neglected.  Ihe  assuniptiaas  made  are 
atpproxlamted  hy  certain  laboratory  plasmas  and  hence  the  results  may  be 
of  practical  iaqportanee. 

Tbr  the  sero  field  case  the  modes  with  azimuthal  variation  are 
found  under  certain  eondlticns  to  edilblt  backward  wave  characteristics 
and  previously  published  results  of  a  quaci-static  analysis  are  modified. 
In  the  ease  of  large  magnetic  field,  the  occurrence  of  unique  vaveguldlng 


features  is  discussed. 
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I,  INOBODUCTIOII 


Ibis  study  is  concerned  with  the  high  frequency  inodes  of  guided 
electronagnetlc  propagation  edong  a  circular  plasma  column  of  Infinite 
length  and  finite  radius,  h,  as  shown  In  Hg.  1.  The  plasma  density  Is 
homogeneous  for  r  <  h  and  decreases  dlscontlnuously  to  zero  at  r  «  b 
where  a  free  spcuie  region  begins.  There  may  also  be  present  a  d.c. 
axial  magnetic  field, 

Specifically  the  cases  of  o  and  oo  are  investigated  In  detail, 
modes  with  azlmuthally  varying  fields  being  Included  In  the  investigation. 
Although  It  has  long  been  recognized  that  Maxwell's  equations  In  these 
two  extremes  of  magnetic  field  are  tractable,  previous  Investigation  Into 
the  mode  structure  of  an  unenclosed  column  has  been  Incomplete. 

Thlvelplece  and  Oould^  have  presented  an  exact  treatment  of  the 

2 

azlmuthally  synmetric  node,  and  TM.velpiece  has  presented  a  quasl-statlc 
treatment  for  the  mode  of  one  angular  variation  for  the  unenclosed  column 
In  zero  magnetic  field.  The  existence  of  backward  waves  and  cutoff  was 
Indicated  for  the  azlmutbally  varying  mode,  but  the  usefulness  of  these 
conclusions  Is  limited  by  the  nature  of  the  quasl-statlc  assumptions. 

Idnhart^  has  examined  the  dispersion  curves  for  the  circularly 
symsKtrlc  modes  In  the  limit  of  large  magnetic  field. 

This  paper  begins  with  a  tutorial  derivation  of  some  basic  equations 
governing  high  frequency  wove  propagation  In  a  homogeneous  plasma  for  the 
general  case  of  arbitrary  magnetic  field.  Subsequent  sections  apply  these 
equations  to  the  derivation  of  dispersion  curves  for  the  two  extreme 
eases  under  Investigation. 


-  3  - 

n.  BASIC  BcajAnoHS 

A.  Capacitivity  Tensor 

A  plaama  may  be  represented  by  e  tensor  capacitivity  whose  substitu¬ 
tion  Into  Maxwell's  equations  yields  the  field  defining  equations.  The 
use  of  a  tensor  of  relatively  simple  form  Involves  the  adoption  of  an 
Idealised  plasma  model.  The  plasma  is  assumed  to  consist  of  a  cold 
coUlslonless  electron  gas  of  density,  n,  Imbedded  In  a  stationary  back¬ 
ground  of  neutralising  Ions  In  the  presence  of  an  axial  magnetic  field, 
Then  with  linearisation  appropriate  to  small  signals  the  motion  of 
the  electrons  Is  given  slaqply  as  the  response  to  a  Lorents  force 


mv  ■  -elE  +  V  X  B^] 


(1) 


where  m  Is  the  electron  mass 

e  Is  the  magnitude  of  electron  charge 
V  is  the  electron  velocity. 

By  defining  a  x>olarlsation  current  density  £  >  n^  and  assuming  time 
dependence  as  e^"^,  the  polarisation  P  may  be  easily  determined  from  (l). 
Then,  setting  0  •  a  *  E  -  9^  +  P  defines  the  tensor  capacitivity  «  . 
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where  Tg  are  subscripts  denoting  arbitrary  ortbogonal 
coordinates  transverse  to  the  axis 


■(^) 


,2\l/2 


Is  the  radian  electron  plasma  frequency 


sBq 

<a^  ■  is  the  radian  electron  cyclotron  frequency. 


B.  Ihe  Fields  of  a  Magneto  Pleana 

Conalder  Maxwell's  equations  In  general  orthogonal  cylindrical  coor- 

2  2  2  2  2  2 

dlnates  with  line  element  ds  ■  h^  dq^  hg  dq^  +  dz  in  a  hosiogeneous 
medium  whose  capacitivity  has  the  form  of  (2)  and  whose  permeability  equals 
the  free  space  value,  Uq*  equations  may  be  written 

(7i  -  W.)  *  (St  *  V.>  *  <’> 

V.)  ■  (fe ’Sr *  •j'A)  W 

Where  T  denotes  the  transverse  part  and  axial  variation 
as  e'^^  is  assumed. 

Equating  the  transverse  parts  In  (?)  and  {k)  yields  the  dependency 
equations  which  esqjnress  the  transverse  field  components  In  terms  of  the 
axial  components 

St  ■  -*%  4'^  ’^T  “  V.)-  h  ‘5) 

(6) 

Where  A  -  Iod^Uq  J  1 


(7) 
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It  renaliiB  to  find  wave  type  equations  governing  the  axial  field 
cGsqponents.  We  proceed  as  is  usual  hy  taking  the  curl  of  (3)  and  substi> 
tuting  fxm  (4)  which  gives 

7x  W  B  -  (E^+  +  (7j  (7y  -JP^)  •  (Sj, 


-  ®^H0  («T  *  ^  +  *3Vs> 


(8) 


Bquating  the  axial  cosponents  of  (8)  gives 

^  \  *  “Vj*.  -  •  5t  (9) 

In  a  similar  manner  hy  taking  the  curl  of  (4)  one  may  show  that 

vj  H.  -  -  -ja>  Wi » •  v  -4  (10) 


The  ejgpressions  on  the  right  of  (9)  and  (lO)  must  new  he  ejqpressed  in 
terms  of  the  axial  components.  To  get  an  eaqpression  for  V*  we  consider 
■Ince  this  is  knew  to  equal  V*  D  >  0.  In 

the  notation  of  general  curvilinear  coordinates 


5r>  -  ^[1^  <VA  *  JVA)  *  ^  (-JVA  *  VA)] 


so  that 


(n) 
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In  a  ainllAr  aunner  ve  nay  evaluate 

fit  c 

(if (“) 

Then  eubatltuting  (11)  into  (9)  and  (12)  into  (lO)  produces  the 
desired  wave  type  equations 


^  (“Vj-  i? if  ®. 

(13) 
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V®  H.  *  («%  .  f?)  .  -Jto  B, 

(lU) 

In  circular  cylindrical  coordinates  with  azimuthal  variation  as 

-1n^ 

e"**  ^  the  dependency  equations  (5)  and  (6)  hecoine 
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and  the  coupled  wave  equations  in  the  cudal  ccaponents  become 
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(15) 


(16) 


(17) 


(18) 


Equations  equivalent  to  (I5),  (16),  (I7)  and  (18)  have  previously 
appeared  in  the  literature/  They  are  in  general  applicable  to  any 
medium  with  tensor  capacitivity  of  the  form  shown  in  (2)  and  by  a  scalar 
permeability.  In  particular,  we  will  apply  them  to  a  plasma  in  two 
limits  of  magnetic  field  where  degenerate  forms  of  the  capacitivity 
tensor  hold  and  the  wave  equations  become  uncoupled. 

These  field  equations  are  also  of  course  applicable  to  the  very 
degenerate  case  of  free  space,  and  we  next  proceed  to  use  them  in  finding 
the  fields  of  the  free  space  region  surrounding  the  plasma  column. 


C.  Solutions  in  the  Free  Space  Region 

In  free  space  the  use  of  equations  (I5),  (16),  (I7)  and  (I8)  entails 

setting  Sgi*  0  and  e^=  6^=  e^.  TSie  wave  equations  \mcouple  and  reduce 

0)^  2 

to  Bessel  equations  with  eigennmber  -  6  .  Hie  appropriate  solutions 

c 

for  bound  modes  are  hyperbolic  Bessel  functions  K^('^)  where 

q  is  real  indicating  that  P  >  ^  • 

The  axial  field  cooiponents  are  then  given  by 

'zo  -  *  "n'f  >  .  B  K„(f ) 

where  A  ar.d  B  are  constants  to  be  determined  by  conditions  at  the 
Interface  r  a  b,  and  by  the  strength  of  excitation. 

The  transverse  fleldsmay  be  found  from  (l^)  and  (16),  the  eiziimithal 


canponents  being 


f, 

V 

j 


-  8  - 

V  -  ^  vf  >  ^  B  #  Vf M  (21) 

L  ^  q  r 

V  [*  ^  ®  ^  ^  Vf  )  (^) 

In  the  next  section  ve  obtain  the  plaana  fields  and  the  associated 
dispersion  relationship  for  a  plasaa  column  in  aero  magnetic  field. 


m  ZERO  MAOimC  FIEID 


As  mentioned  in  the  introduction,  the  ease  of  no  azimuthal  variation 

1  2 

(n  ■  O)  has  been  treated  previously  ’  ,  and  therefore  this  Investigation  is 
confined  to  the  higher  order  nodes  (n  >  1). 


A.  The  Characteristic  Equation 


(a)  Derivation 


With  ■  0  the  capacitivity  tensor  of  equation  (2)  degenerates  into 
o  o 


the  scalar  Cq(1 


),  and  the  wave  equations,  (17)  and  (18),  reduce  to 


Bessel  equations  with  eigennumhers  ~2  ^  ^  ~  B  •  plasma  region 

solutions  are  accordingly 

p2.  ( 

c®  c®  c® 

Note  that  the  factor  is  always  positive  so  that  the  radial 

c= 

eigennumber  in  the  plasma  region,  p,  is  always  greater  than  that  in  the 
free  space  region,  q.  Equation  (24)  is  a  hyperbola  on  the  p,  q  plane 

OCnb 

which  intercepts  the  q  ■  0  axis  at  p  ■  . 

Transverse  field  components  may  be  determined  from  equaticms  (1^) 


and  (16);  the  azimuthal  components  are 


B^i  -  C 


p2r  Vb^  °  p2 
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]jnposlng  the  condition  that  the  tangential  and  z)  field  canqpoaenta 
be  matched  at  the  boundary  r  ■  b  gives  a  set  of  four  homogeneous  equations 
In  k,  B,  C  and  D.  For  nontrivial  solutions  the  determinant  of  the 
coefficient  matrix  must  be  set  equal  to  zero,  which  results  as  usual  In 
the  characteristic  equation.  It  proves  convenient  to  formulate  the 
characteristic  equation  in  terms  of  the  foUowlDg  functions. 


XX  (x)  ' 
n'  ' 


x„  ,(x) 

’ 


i;(x) 


Vi(^> 


(27) 


nie  characteristic  equation  is  then 

2 

[li(p)  +  <(q)]  t(l-  ^  )  i;(p)  +  K;(q)]  +  ti;(p)  +  K;(q)l 

2 

i>)  +  <(q)  1-0  (28). 

Ihe  higher  order  modes  (n  >  1)  governed  by  this  eqiiation  are  hybrid 
possessing  axial  cooponents  of  both  E  and  H.  Ibe  existence  of  ZB  or  IN 
modes  woxxld  imply  n(p  -  q  )  ■  0  which  contradicts  (24). 

The  hyperbolic  Bessel  functions  of  first  and  second  kind  which 
appear  in  eqviatlon  (28)  are  monotonlc,  indicating  the  existence  of  a 
single  mode  for  each  value  of  n. 


(b)  Non-Existence  of  Solutions  for  co  > 

It  may  be  shown  that  no  solutions  exist  for  (28)  when  a>  >  . 

This  Implies  that  for  the  existence  of  guided  waves  there  Is  a  maximum 
frequency  limit  or  equivalently  the  constraint  that  the  capacitivity  ' 
must  be  negative  and  greater  in  magnitude  than  c^. 

The  characteristic  equation  (28)  may  be  written  in  the  form 


2  Kp(q)  K;(q)  +  K^(q)  i;(p)  K;(q)  i^(p) 
2  r(p)  i^(p)  +  kJ;(«i)  i;(p)  +  K;(q)  i*(p) 


Denonstration  of  the  fact  that  a>  <  ^  is  thus  eg.u 
K^(4)  ^  Op)  ^n^P)  PPP  ^  expandad  aa 

.=  Vl(’)  W’)  .  Vl(r) 

-  .  .  _  -  >  "A 


thus  agplvalent  to  shewing  that 


Hote  the  following  Integral  es^esslons  for  the  products  of 

hyperbolic  Bessel  functions a 

00 

K  (q)  K^(q)  ■  2  J  K^^^(2q  cosh  t)  cosh  (|i-k)  t  dt 
0 

*/2 

I  (P)  Iyj(p)  ■  f  J  ®  *** 


From  (51), 


0 

00 

and  K^(q)  -  2  J  Kg^^CBq  cosh  t)  dt  (3 

0 

Now  cosh  2t  >  1,  and  q  is  real  so  that  the  integrands  In  (55)  and  (54)  are 
always  positive.  Ifaus  ve  see 

Vl''*)  2.  4(1)  <= 


12 


In  a  slBilar  Banner  uaing  (32)«  it  aay  be  bbown  that 

>  Vl(P)  Wp)  (5») 

Then  (3?)  anl  (36)  together  with  the  fact  that  p  >  q.  establiah  the 
required  inequality  (30)« 


(c)  Cutoff  and  Resonance 

nie  characteristic  equation  (28)  Is  represented  by  the  solid  lines 

In  Fig.  2  for  n  >  1  and  in  Fig.  3  for  n  ■  2  for  ▼arlous  values  of  the 

a£ 

peraaeter  ^  .  The  features  of  the  curves  for  n  •  2  nay  be  taken  as 

aji 

representative  for  modes  of  larger  n. 

Note  that  the  curves  for  n  >  1  begin  at  the  origin  of  the  p,  q 
plane  whereas  for  n  ■  2  the  curves  intersect  the  q  ■  0  axis  at  Beam  non¬ 
zero  value  of  p.  Also,  both  sets  of  curves  intersect  the  line  p  ■  q  for 

4 

sbbU  values  of  .  Ibese  features  are  of  physical  significance  and  it 

ODF 

is  of  value  to  trace  their  occiurence  mathematically. 

nie  q  «  0  axis  represents  a  cutoff  for  bound  modes;  from  (I9)  it 
may  be  seen  that  at  q>0,  which  is  the  transition  point 

between  slow  and  fast  waves.  The  line  q  ■  p  is  a  resonance  since  from 
(2lf)  it  implies  3  •*  00 . 

At  cutoff  (q  ■  0),  note  the  following  limits  of  the  functions  of  q 
appearing  in  (28) 


(37) 

q  -*  0  ^ 

11m  lC(q)  ■  -In  q,  n  ■  1 

(58) 

1 

0 

I 

o' 

^llm^  I^(q)  -  ,  n  >  1 

(39) 

Figure  2  Cbaracterlstic  Bqvstlon  for  the  Flasna  Column 
In  Zero  Kegnetle  Field,  n  >  1 


equation  (28) 
equation  (2k) 


liAiitow  mi  rtf 


-  15  - 

Beeauee  of  the  distinct  cases  In  (58)  and  (59)  the  n  ■  1  mode  must  he 
treated  separately  from  n  >  1. 

Fbr  the  mode  of  one  angular  variation,  the  characteristic  equation 
in  the  cutoff  limit  hecones 


where  the  finite  nature  of  Iq(p)  has  been  considered.  Equation  (1»0) 
can  be  satisfied  for  ncmzero  a>  only  when  l~(p)  **  co  which  Ingplles 
p  a  0.  Ibis  restilt  will  be  seen  to  lead  to  a  zero  cutoff  frequency  for  the 
n  ■  1  mode. 

For  the  modes  with  n  >  1,  the  characteristic  equation  at  cutoff  for 
nonzero  <0  becomes 


This  can  be  satisfied  only  when 


i;(p) 


.  2)(n.l) 


(U2) 


Ihe  value  of  p  defined  by  (42)  is  the  intercept  on  the  q  ■  0  axis 

which  starts  at  p  ■  0  for  the  marginal  case  ^  »  2  and  rises  to  higher 

or 

values  of  p  with  decreasing  a>.  The  nonzero  value  of  the  Intercept  will 
be  seen  to  imply  a  nonzero  cutoff  frequency  for  modes  with  n  >  1. 


I 
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Next  consider  the  resonance  limit  q  ■  p.  Ihe  characteristic 
equation  (28)  hecomes 

<Vl  *  Vl  V  f  1 

"  »n-l  *»  *  Vl  ( t"-  ^  1  Vl  *  Vl  V  -  0  («) 

Where  the  argument  of  both  the  I  and  K  functions  is  p. 

The  identity 

Vi(p)  -  5  (■*> 

may  be  used  to  reduce  (^3)  to  the  form 

This  relationship  is  sketched  in  figure  k. 

It  will  be  noted  that  for  small  equation  (43)  does  not  have  a 

^  Op 

solution.  However^  for  values  of  o  approaching  there  are  two 

o  ^ 

values  of  p  which  satisfy  (43)  for  a  given -  .  It  will  subsequently 

“p 

be  shown  that  these  soluticms  indicate  the  presence  of  bcuikward  waves. 


B.  The  Dispersion  Curves 


Substituting  for  p  and  q  from  (24)  and  (I9)  into  the  character- 

f  A 

Istlc  equation  (28)  results  in  a  dispersion  relationship  between  ~  and 

Opb 


0b  with 


as  a  parameter. 


[i;(i) .  <(«)]  (if)  1  ♦  'it*)] 

.  [i;(p) .  isjci)]  (?)  i>) .  <(«)]  -  0 

, .  [(pb)=.  p  -  [((«.)=-  (f  )=♦ 


>  (w) 


y 


p 


Characteristic  Eg^iatlon  Resonance  Condltlan 
Bjostlon  (If^) 
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Points  satisfying  this  equation  may  be  obtained  by  finding  the 
intersection  of  a  hyperbola 


2  2 
p  -  q 


(2k) 


with  p"  ots  of  the  characteristic  equation  for  various  values  of 


A  point  of  intersection  gives  q,  p, 


0) 


and 


A 

'  (D 


Which  yields 


uib 

c 


(^)(^)  »d  Sb  .  . 


From  the  hyperbolas  for 


2^  and 


2Eb  _  _ 


7^  plotted  as  dashed 


c  '  c 

lines  in  figure  2,  it  nay  be  seen  that  the  dispersion  curves  for  the  mode 

of  one  angular  variation  (n  >  1)  passes  through  the  origin,  ~  ■  0, 

c 


and  asynqptotlcally  approaches 
of 


1  ^ 


as  ^  CO .  For  small  values 


tupb 


c  ^  c 

the  dispersion  ciurves  will  be  triple  valued  in  3  for  a  single 


value  of  oj  as  indicated  by  the  triple  intersection  of  the  hyperbola  for 
m  .25  with  the  characteristic  equation  curve  for  ^  *  2.};  this 
implies  the  existence  of  a  region  of  negative  slope  or  backward  wave 
propagation. 


hyperbolas  for  =  .5  and  ■1.5  are  represented  by  dashed 


ob  ^  COnb 

=  .5  and  - 
c  ^  c 

lines  in  Pig.  5.  Ihe  nature  of  their  intersection  with  the  plots  of  the 


characteristic  equation  for  n  ■  2  shew  that  for  n  >  1  the  dispersion 
curves  will  begin  on  the  line  —  ■  ^b  at  some  point  other  than  the  origin, 
nils  cutoff  frequency,  is  determined  from  (k2)  as 


cutoff 


a>_ 


2(n-l)  .  1, 


1/S 


m 


uib  1  ain^ 

nie  dispersion  curves  asynqptotically  approach  ~  ■ 


values  of 


For  small 


,  the  dispersion  curves  will  be  double  valued  in  p  again 


implying  the  existence  of  backward  wave  propagation. 
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Dispersion  curves  for  n  ■  0  are  reproduced  from  previously 
published  data^  in  Fig.  ^  for  purposes  of  cooparlson.  Dispersion  cvirves 
for  n  ■  1  are  plotted  in  Fig.  6  and  for  n  a  2  in  Fig.  7  for  revealing 
values  of  the  parsmeter  .  As  indicated  above,  an  important  difference 

between  the  curves  for  the  n  a  o  mode  and  the  modes  of  higher  order  is 
the  appearance  of  backward  waves  in  the  azlnmthally  varying  modes.  The 
portions  of  the  cmnres  exhibiting  this  behaviour  are  enlarged  in  the  lower 
right  comers  of  Pig.  6  and  7*  TTxe  phenomena  of  backward  waves  must  be 
associated  with  either  the  properties  of  the  plasma  as  a  negative 
dispersive  permittivity  or  with  the  hybrid  nature  of  the  higher  order 

7 

modes  or  both. 

The  disappearance  of  the  backward  wave  characteristics  for  large 

OJnb 

values  of  should  be  noted.  Increasing  this  parameter  Increases  p 

relative  to  q  [see  (24)]  and  this  has  the  effect  of  decreasing  field 

penetration  into  the  body  of  the  plasma  and  of  concentrating  power  flow 

in  the  free  space  region  as  may  be  seen  from  consideration  of  the  nature 

of  the  hyperbolic  Bessel  fimctions. 

All  three  sets  of  dispersion  curves,  Fig.  5,  6  and  7,  show  resonance 

at  o)  a  ^  ,  i.e.,  p  increases  without  limit  as  this  frequency  is 

\JZ 

approe^hed.  Increasing  0  Implies  Increasing  p  and  q  and  consequently 
increasing  confinement  of  field  energy  to  the  region  of  the  interface. 
I\irthermore,  group  velocity  in/ap  rapidly  approaches  zero  at  the 
resonance  so  that  energy  does  not  propagate. 

As  0  Increases,  the  phase  velocity,  (o/0,  of  course  decreases,  and 
when  its  value  approaches  the  electron  thermal  velocity,  the  assumption 
of  negligible  thermal  pressure  made  in  section  II  becomes  tintenable.  An 


so  - 


0  12  3  4  5  6 


/3b 

Figure  ^  Dispersion  Curves  for  Plasaa  Cblum 
with  Sero  Magnetic  Field,  n  >  0 
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Ueperelon  Curves  for  Plasma  Column 
with  Zero  Nsgnetlc  field,  n  ■  1 


i 

i 

4 

I 

i 

i 


I 

5 


I 


i 


Rigorous  Analysis 
Quasl-Statie  Analysis 


! 

I 


-  23  - 


analysis  that  Includes  the  effect  of  thermal  pressure  would  indicate 

that  as  phase  velocity  approaches  electron  thermal  velocity,  damping 
8  Q 

set  in  ,  the  field  energy  being  converted  to  electron  kinetic  energy. 


C.  Comparison  with  Quasi-Static  Analysis 


It  is  of  Interest  to  ccmpare  the  modal  properties  presented  above 

with  the  predictions  of  approximate  analysis  using  the  quasi-static 

2 

assumption  as  carried  out  by  Trlvelplece  for  the  n  >  1  mode. 

nie  basic  quasi-static  assumption  is  -VV,  V  being  a  scalar 

potential.  For  the  unmagnetized  plasma  column  this  assumption  implies 
2  a? 

not  only  that  p  »  with  non -dispersive  media  but  also  that 

2  “I 

e  »3' 

Ihe  quasi-static  assumption  leads  to  the  following  dispersion 


equation 


2 

(U 


K;(pb)  +  K^(pb) 

I"  Ob)  +  I^Ob) 


Equation  (U8)  for  the  n  ■  1  mode  is  represented  by  the  dashed  lines  in 

Pig.  6.  Hiey  have  intercepts  with  the  p  »  0  axis  at  o)  »  ^  and 

V2 

corresponding  nonzero  intercepts  with  the  line  p  *  ^  .  Uie  region  of 

these  intercepts  (P<~)  is  not  with  the  province  of  quasi-static 

analysis  and  their  existence  was  correctly  viewed  with  suspicion  by 

Trlvelplece.  His  curves  also  indicate  backward  waves  for  small  pb,  a 

feature  which  is  retained  in  the  exact  analysis  only  for  small  values  of 
oupb 

c  * 


(U8) 


the  parameter 
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In  the  next  section  ve  consider  the  plasma  column  In  large  axial 
magnetic  field.  Modal  properties  In  that  case  will  he  seen  to  differ 
significantly  from  those  presented  above  for  the  opposite  extreme  of 
the  xinmagnetlzed  column. 
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rv.  THE  PLASMA  COLUMN  IN  LARGE  AXIAL  MAGNETIC  FIELD 


Large  Eixlal  magnetic  field  is  meant  to  Imply  an  electron  cyclotron 
frequency  greatly  in  excess  of  either  the  electron  plasma  frequency  or 
the  signal  frequency,  llils  condition  is  not  tinrealistic  for  laboratory 
plasmas.  Often  a  laboratory  plasma  has  an  axial  magnetic  field  variable 
up  to  20  kilogauss,  which  corresponds  to  a  cyclotron  frequency, 
f^  =  56  gc./sec.,  while  a  typical  value  of  electron  density  might  be 
10"  cm.”^  which  corresponds  to  a  plasma  frequency,  f^  «  2.81*  gc./sec.; 
it  will  be  shown  that  for  guided  propagation  the  signal  frequency  is 
smaller  than  f^.  The  large  axial  magnetic  field  will  not  Interact  with 
electron  motion  parallel  to  it  but  will  effectively  inhibit  any  transverse 
motion. 

A.  The  Dispersion  Curves 


In  the  limit  of  large  B^,  the  capacitivity  tensor  of  equation  (2) 


becomes  uniaxial 


e  -•  €, 


1 

0 

0 


0 

1 

0 


0 

0 

.A 
^  2 
cu 


m 


The  wave  equation  (17.)  and  (18)  uncouple,  and  the^equation  in  E^  trans¬ 
forms  to  Bessel's  equation  with  eigennumber  (1-^)  while  the 

“  °  a,2  g 

equation  in  beccxnes  Bessel's  equation  with  eigennumber  (^  -  P  )• 

Consequently,  H  is  unaffected  by  the  plasma  and  no  guided  TE  modes 
2 


are  possible 


-  26  - 


Hie  axial  conqponent  of  E  will  have  the  fora 


E  -  F  J  (^) 
zl  n'  b  ' 


4  -  4-1) 0^-4.=)  =  4-1),= 


£J3“ 


(50) 

(51) 


For  bound  modes  p_  is  real  imposing  the  constraint  o)  <  a>  . 

lU  P 

Hie  modes  are  H4  having  field  components  E  ,  E  ,  Hj  for  the 

z’  r'  jD 

n  ■  0  case  and  all  field  components  except  H  for  n  >  1.  Hie  charae- 
teristic  equation  may  be  found  as  outlined  in  the  previous  section  and  is 


J  (P_) 
_  n  ni 

^m  J  (p  ) 


<(,) 
1  w 


(52) 


Hiis  equation  has  previously  appeared  in  the  literature.^® 

Hie  left  side  of  (52)  is  oscillatory  and  there  will  consequently  be 
an  infinity  of  solutions  for  a  given  n.  Hie  modes  are  then  doubly 
infinite  in  number  and  a  given  mode  will  be  designated  by  a  double  sub¬ 


script  as  Ht 


nm 


Hie  dispersion  relationship  is  gotten  by  substituting  in  the 

Q)b 

characteristic  equation  for  p  and  q  in  terms  of  —  ,  Sb,  and 

m  0 

oi^b 

as  determined  by  equctions  (I9)  and  (51).  Hie  dispersion  relation 

is  plotted  for  several  modes  in  Fig.  8  with  s  10. 

c 

Hie  curves  have  similar  shapes  for  different  modes;  they  all  begin 
at  the  origin  and  asymptotically  approach  the  line  o  ■  (Pp  for  large  6. 
As  in  the  zero  field  case  the  asymptotic  approach  is  accompanied  hy 
resonant  excitation  of  electron  oscillations.  Along  the  asymptote,  q 


Increases  without  bound,  but  the  plasma  elgennumber  p  remains 

m 
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essentially  constant,  and  the  fields  are  concentrated  In  the  body  of  the 

plasma  rather  than  at  the  Interface. 

Tig.  8  also  demonstrates  that  the  waves  heccane  slower  (l.e.,  phase 

and  group  velocity  decrease)  as  the  order  of  the  mode  Increases,  nie 

initial  slope  of  the  dispersion  curve  may  be  shown  to  equal 
I-  (n_b  _2  -I  1/2 

1  +  (-^)  z  .  for  the  m  mode,  where  z  .  is  the  mth 

L  '  c  '  n-l,i^  nm  '  n-l,m 

zero  of  J  ^  exclusive  of  the  zero  at  the  origin  except  for  n  >  0. 
n-i 

nie  initial  slope  of  the  mode  is  unity.  Trlvelpiece  and  Oould^ 

found  ctirves  similar  to  Fig.  8  for  a  plasma  filling  a  circular  conducting 
pipe.  In  that  case  the  initial  slopes  of  the  curves  were  found  to  be 

so  that  the  waves  are  seen  to  be  faster  for  the  open 

column. 


B.  A  Unique  Vtevegulding  Feature 

There  is  a  rather  vuiique  feature  of  the  equations  for  this  case  of 
large  magnetic  field  which  has  interesting  physical  implications. 

Equation  (^)  expresses  In  terms  of  q  indei>endently  of  column 
radius,  b,  and  the  characteristic  equation  (^)  also  defines  a  relation¬ 
ship  between  the  radial  eigennumbers  which  is  independent  of  b.  Thus 
either  q  or  p  may  be  determined  as  a  function  of  the  ratio  ~  for 

®  ojp 

arbitrary  column  radius.  Curves  of  q  vs.  ^  are  plotted  in  Pig.  9 

®P 

for  several  modes.  That  radial  eigenvalues  are  Independent  of  the  trans¬ 
verse  of  the  column  may  be  associated  with  the  fact  that  the 

portion  of  the  capacitivity  tensor  pertaining  to  the  coordinates  transverse 
to  the  magnetic  field  is  a  scalar  equal  to  c^. 
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The  effect  on  wavegulding  properties  of  the  colutm  nay  be  seen  by 
considering  equation  (19)  which  we  rewrite  as 


3" 


+  a! 


(53) 


For  a  fixed  and  o)  (and  therefore  a  fixed  q)  as  the  column  radius  is 
increased  phase  velocity  approaches  the  speed  of  light  in  vacuo  (i.e., 

3  ~  •  Pig*  10  depicts  the  dependence  of  the  dispersion  curve  for  the 

mode  on  coltimn  radius  b. 

Bie  tendency  of  the  waves  to  become  faster  with  increasing  radius  is 
opposite  to  the  behaviour  of  an  open  dielectric  rod  waveguide.  In  the  latter 
case,  if  the  radius  of  the  dielectric  rod  greatly  exceed  the  wavelength  in 
the  dielectric,  phase  velocity  is  slowed  to  the  speed  of  li£dit  in  the 


dielectric. 


.  On  the  other  hand,  if  the  radius  of  the  plasma  column 


is  much  larger  than  the  wavelength,  we  approach  the  situation  of  having  all 
the  field  energy  stored  in  field  components  transverse  to  the  direction  of 
propagation  (except  at  resonance)  so  that  the  effective  capacitivity  is 
Cq  and  phase  velocity  approaches  c. 

To  calculate  stored  energy  one  uses  the  expression  for  average  energy 
density  appropriate  to  dispersive  medic 

3(aj«  ) 


_  1  f  3(aj«  )  \ 

U  =  ^  Real  B  +  •  e) 


(54) 


Then  by  substituting  into  (54)  the  appropriate  field  expressions  for  the 
mode  under  consideration  which  may  be  obtained  from  equations  (50),  (20), 
(15),  and  (16),  and  by  integrating  U  over  the  cross-section  one  may  arrive 
at  expressions  for  stored  energy.  Specifically,  the  ratio  of  energy  in 
the  transverse  fields  to  that  in  the  axial  field  for  the  TML^  mode  is 


I 


where  the  argument  of  the  J's  is  p 

m 

and  the  argument  of  the  K's  is  q 

Equation  makes  clear  the  increasing  storage  of  energy  in  the  transverse 
fields  with  increasing  h. 

As  the  signal  frequency  approaches  the  plasma  frequency  the  augiaented 
axial  oscillations  of  the  electrons  near  resonance  is  accompanied  by 
increasing  storage  of  energy  in  the  axial  field  and  the  wave  is  slowed. 


i. 


i 

i 

i 
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V.  CONCLUSIQWS 

The  properties  of  guided  electromagnetic  propagation  along  an  vinen- 
closed  circular  cylinder  of  plasma  have  been  shown  to  contrast  markedly 
for  the  two  extremes  of  axial  magnetic  field. 

For  the  case  of  zero  magnetic  field  a  single  bound  mode  was  found 
for  each  value  of  the  Integer  n  where  the  azimuthal  variation  of  the 
modes  is  as  e  .  Uie  azimuthally  varying  modes  (n  >  1)  are  hybrid. 
Hie  mode  of  one  angular  variation  has  no  cutoff  and  its  modal  properties 
predicted  by  quasi-static  analysis  are  significantly  modified  by  the 
results  of  a  more  rigorous  treatment.  Modes  of  n  >  1  can  exist  only 
over  a  narrow  range  of  frequency  having  a  nonzero  value  of  cutoff 
flrequency.  All  modes  experience  a  resonance  at  ui  «  ^  and  are  nonexist- 

\/l 

ent  at  higher  frequencies.  The  plasma  fields  vary  radially  as  hyperbolic 
Bessel  functions  of  the  first  kind,  I^.  Perhaps  the  most  significant 
feature  is  the  appearance  of  backward  waves  for  the  azimuthally  varying 

(jDpb 

modes  for  small  values  of  the  parameter  . 

With  large  axial  magnetic  field  there  is  an  infinity  of  TM  modes 

for  each  value  of  n.  All  modes  begin  at  zero  frequency  and  experience 

a  resonance  at  cu  >■  (u  being  nonexistent  for  w  >  w  .  The  fields  in 

P  P 

the  plasma  region  vary  radially  as  Bessel  functions  of  real  argument, 

J  .  Hie  radial  elgennumbers  are  independent  of  colvmin  thickness,  and 
n 

phase  velocity  approaches  the  speed  of  light  in  vacuo  if  the  column  radius 
is  large  canpared  to  the  wavelength. 
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